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HOLOMORPHIC VECTOR FIELDS AND QUADRATIC 
DIFFERENTIALS ON PLANAR TRIANGULAR MESHES 

WAI YEUNG LAM AND ULRICH PINKALL 


Abstract. Given a triangulated region in the complex plane, a discrete vector 
field Y assigns a vector Yi G C to every vertex. We call such a vector field 
holomorphic if it defines an infinitesimal deformation of the triangulation that 
preserves length cross ratios. We show that each holomorphic vector field 
can be constructed based on a discrete harmonic function in the sense of the 
cotan Laplacian. Moreover, to each holomorphic vector field we associate in a 
Mobius invariant fashion a certain holomorphic quadratic differential. Here a 
quadratic differential is defined as an object that assigns a purely imaginary 
number to each interior edge. Then we derive a Weierstrass representation 
formula, which shows how a holomorphic quadratic differential can be used to 
construct a discrete minimal surface with prescribed Gaufi map and prescribed 
Hopf differential. 


1. Introduction 


Consider an open subset U in the complex plane C = R 2 with coordinates 
z = x + iy together with a holomorphic vector field 

Y = f —. 

J dx 

Here Y is a real vector field. It assigns to each pel 2 the vector f(p) £ C = K 2 . 
We do not consider objects like ^ which are sections of the complexified tangent 
bundle (TR 2 ) C . 

Note / : U — > C is a holomorphic function, i.e. 


0 = f. 


1 (d_l 

2 \dx 



Let t H > g t denote the local flow of Y (defined for small t on open subsets of U 
with compact closure in U). Then the euclidean metric pulled back under g t is 
conformally equivalently to the original metric: 

gU,)=e 2u (,) 


for some real-valued function u. The infinitesimal change in scale u is given by 

u = ^div Y = Re (/ z ). 

Note that u is a harmonic function: 


Uzz = 0. 

On the other hand, differentiating u twice with respect to z yields one half the third 
derivative of /: 

'il'ZZ = 2 JzZZ‘ 
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It is well-known that the vector field Y corresponds to an infinitesimal Mobius 
transformation of the extended complex plane C if and only if / is a quadratic 
polynomial. In this sense f zzz measures the infinitesimal “change in Mobius struc¬ 
ture” under Y (Mobius structures are sometimes also called “complex projective 
structures” 0). Moreover, the holomorphic quadratic differential 

Q ■= fzzz dz 2 

is invariant under Mobius transformations d>. This is equivalent to saying that q is 
unchanged under a change of variable $(;:) = w = £ + ig whenever >I> is a Mobius 
transformation. This is easy to see if = az + b is an affine transformation. In 
this case 


dw = adz 


and therefore 


with 


d Id 
dw a dz 


Y = 



/ = «/■ 


Thus we indeed have 


f dvr — f d 2 T 

J WWW UjUJ — J zzz • 


A similar argument applies to d>(z) = f and therefore to all Mobius transformations. 

For realizations from an open subset U of the Riemann sphere CP 1 the vanishing 
of the Schwarzian derivative characterizes Mobius transformations. The quadratic 
differential q plays a similar role for vector fields. We call q the Mobius derivative 
of y. 

An important geometric context where holomorphic quadratic differentials arise 
comes from the theory of minimal surfaces: Given a simply connected Riemann 
surface M together with a holomorphic immersion g : M — > S 2 C R 3 and a holo¬ 
morphic quadratic differential q on M , there is a minimal surface F : M —> R 3 
(unique up to translations) whose Gaufi map is g and whose second fundamental 
form is Req. 

In this paper we will provide a discrete version for all details of the above story. 
Instead of smooth surfaces we will work with triangulated surfaces of arbitrary 
combinatorics. The notion of conformality will be that of conformal equivalence as 
explained in (3j. Holomorphic vector fields will be defined as infinitesimal conformal 
deformations. 

There is also a completely parallel discrete story where conformal equivalence of 
planar triangulations is replaced by preserving intersection angles of circumcircles. 
To some extent we also tell this parallel story that belongs to the world of circle 
patterns. 

The results on planar triangular meshes in this paper are closely related to 
isothermic triangulated surfaces in Euclidean space [8]. 


2. Discrete conformality 

In this section, we review two notions of discrete conformality for planar trian¬ 
gular meshes. We first start with some notations of triangular meshes. 

Definition 2.1. A triangular mesh M is a simplicial complex whose underlying 
topological space is a connected 2-manifold (with boundary). The set of vertices 
(0-cells), edges (1-cells) and triangles (2-cells) are denoted as V, E and F. 
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We denote E int the set of interior edges and V mt the set of interior vertices. 
Without further notice we will assume that all triangular meshes under considera¬ 
tion are oriented. 

Definition 2.2. A realization z : V —> C of a triangular mesh M in the extended 
complex plane assigns to each vertex i £ V a point Zi £ C in such a way that 
for each triangle {ijk} £ F the points corresponding to its three vertices are not 
collinear. 

Given two complex numbers z\ , Z 2 £ C we write 

(zi,z 2 ) ■= Re( 2 i 2 2 )- 

We are looking for suitable definitions of conformal structure of a realization z. 
In particular, we want z to be conformally equivalent to g o z whenever g : C —» C 
is a Mobius transformations. This requirement will certainly be met if we base 
our definitions on complex cross ratios: Given a triangular mesh z : V —> C, we 
associate a complex number to each interior edge {ij} £ Ei nt , namely the cross 
ratio of the corresponding four vertices (See Figure [l]) 

_ (Zj - Zk){Zj - Zi) 

CTz ’ lJ (Zk - Zi)(zi - Zj)' 

Notice that cr Z jj = cr z ji and hence cr 2 : Ei nt —> C is well defined. It is easy to 



Figure 1. Two neighboring and oriented triangles 

see that two realizations differ only by a Mobius transformation if and only if their 
corresponding cross ratios are the same. In order to arrive at a more flexible notion 
of conformality we need to relax the condition that demands the equality of all cross 
ratios. Two natural ways to do this is to only require equality of either the norm or 
alternatively the argument of the cross ratios. This leads to two different notions of 
discrete conformality: conformal equivalence theory png and circle pattern theory 

IIIl- 

Note that for the sake of simplicity of exposition we are ignoring here realizations 
in C where one of the vertices is mapped to infinity. 

2.1. Conformal equivalence. The edge lengths of a triangular mesh realized in 
the complex plane provide a discrete counterpart for the induced Euclidean metric 
in the smooth theory. A notion of conformal equivalence based on edge lengths was 
proposed by Luo [9]. Later Bobenko et al. [3] stated this notion in the following 
form: 

Definition 2.3. Two realizations of a triangular mesh z,w \ V —¥ C are conformally 
equivalent if the norm of the corresponding cross ratios are equal: 

| cr z | = | cr w |, 

i.e. for each interior edge {ij} 

I (Zj - Z k )\\(Zj - Zi) I _ I (Wj - w k )\\(Wj - wj) I 
lOfc - Zi)\\{zi - Zj )I I (w k - Wi)||(wi - Wj)\' 
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This definition can be restated in an equivalent form that closely mirrors the 
notion of conformal equivalence of Riemannian metrics: 

Theorem 2.4. Two realizations of a triangular mesh z,w : V —> C are conformally 
equivalent if and only if there exists u : V —> K. such that 

u i+ u j 

I Wj - wt | = Zj - Zi\. 

Proof. It is easy to see that the existence of u implies conformal equivalence. 
Conversely, for two conformally equivalent realizations z, w, we define a function 
a : E —> K. by 

|Wj - Wi\ = e aij \zj - Zi\. 

Since z,w are conformally equivalent cr satisfies for each interior edge {ij} 

&jk &ki T Cn (7ij = 0. 

For any vertex i and any triangle {ijk} containing it we then define 

gM; gCfci+fy — <?jk 

Note the vertex star of i is a triangulated disk if i is interior, or is a fan if i is 
a boundary vertex. Hence the value Ui defined in this way is independent of the 
chosen triangle. □ 

2.2. Circle patterns. Given a triangular mesh realized in the complex plane we 
consider the circumscribed circles of its triangles. These circles inherit an orienta¬ 
tion from their triangles. The intersection angles of these circles from neighboring 
triangles define a function f) : Ei nt —> [0, 27r) which is related to the argument of 
the corresponding cross ratio via 

(1) = Arg(cr Z) jj). 



Figure 2. The intersection angle of two neighboring circum¬ 
scribed circles 

Based on these angles we obtain another notion of discrete conformality which 
reflects the angle-preserving property that we have in the smooth theory. 

Definition 2.5. Two realizations of a triangular mesh z,w : V —>■ C have the same 
pattern structure if the corresponding intersection angles of neighboring circum¬ 
scribed circles are equal: 

Arg(cr Z) jj) = Arg(cr. u , > i :/ ), 
i.e. for each interior edge {ij} 

Arg fa - ~ zi) = Arg {Wj - W k )[Wj - wi) 

8 {z k - Zi)(zi - Zj) B (w k - Wi)(wi - Wj) ’ 

Just as conformal equivalence was related to scale factors u at vertices, having 
the same pattern structure is related to the existence of certain angular velocities 
a located at vertices: 
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Theorem 2.6. Two realizations of a triangular mesh z,w : V —> C have the same 
pattern structure if and only if there exists a : V —> [0, 2ir ) such that 

w-i — Wi ..“*+“3 Zn — Zi 

_ ± _= e 2 _-_. 

\Wj-Wi\ \Zj-Zi\' 

Proof. The argument is very similar to the one for Theorem |2.4| In particular, 
the existence of the function a easily implies equality of the pattern structures. 
Conversely, assuming identical pattern structures we take any oj : E —> R that 
satisfies 

Wj - m = Zj - Zj 

\Wj-Wi\ \Zj-Zi\' 

For any vertex i and any triangle {ijk} containing it we define cq £ [0, 27 t) such 
that 

e iai _ gjfwfei+wy-wjfc) 

Note the vertex star of i is a triangulated disk if i is interior, or is a fan if i is a 
boundary vertex. Hence having the same pattern structure implies that the value 
on is independent of the chosen triangle. □ 

3. Infinitesimal deformations and linear conformal theory 

We will linearize both of the above notions of discrete conformality by considering 
infinitesimal deformations. This will allow us to relate them to linear discrete 
complex analysis, based on a discrete analogue of the Cauchy Riemann equations 
[3HHD] (See the survey jl2]h 

Definition 3.1. An infinitesimal conformal deformation of a realization z : V —> C 
of a triangular mesh is a map z : V — > C such that there exists u : V —> R satisfying 

Zj - Zj = {.Zj - Zj , Zj - Zj) Uj + Uj 

Zj-Zi \zj - Zi\ 2 2 • 

We call u the scale change at vertices. 

Definition 3.2. An infinitesimal pattern deformation of a realization z : V —> C of 
a triangular mesh is a map z : V —> C such that there exists a : V —>■ R satisfying 

l m ~ £i = fa - Zj,i(zj - zj) = + a 0 

Zj-Zi \zj - Zi \ 2 2 

We call a the angular velocities at vertices. 

Example 3.3. The infinitesimal deformations z := az 2 + bz + c, where a, 6, c £ C 
are constants, are both conformal and pattern deformations since 

——^ = (azt + 6 / 2 ) + (azj + 6 / 2 ). 

Zj 

Infinitesimal conformal deformations and infinitesimal pattern deformations are 
closely related: 

Theorem 3.4. Suppose z : V —► C is a realization of a triangular mesh. Then 
an infinitesimal deformation z : V —> C is conformal if and only if iz is a pattern 
deformation. 

Proof. Notice 

{zj - Zj, Zj - Zj) = (izj - izj,i{zj - zj) 

\Zj-Zi\ 2 \Zj ~ Zi I 2 

and the claim follows from Definition 13.11 and 13.21 □ 
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3.1. Infinitesimal deformations of a triangle. Let z : V C be a realization 
of a triangulated mesh and z an infinitesimal deformation. Up to an infinitesimal 
translation z is completely determined by the infinitesimal scalings and rotations 
that it induces on each edge. These infinitesimal scalings and rotations of edges 
satisfy certain compatibility conditions on each triangle. These conditions involve 
the cotangent coefficients well known from the theory of discrete Laplacians. As we 
will see in section 3.2 for conformal deformations (as well as for pattern deforma¬ 


tions) the infinitesimal scalings and rotations of edges are indeed discrete harmonic 
functions. 

Consider three pairwise distinct points Z\, z 2 , 23 £ C that do not lie on a line. 
In the following i,j,k denotes any cyclic permutation of the indexes 1,2,3. The 
triangle angle at the vertex i is denoted by /3*. We adopt the convention that 
all /3 i,/3 2 ,/3 3 have positive sign if the triangle Zi,Z 2 ,z 3 is positively oriented and 
a negative sign otherwise. Suppose we have an infinitesimal deformation of this 
triangle. Then there exists aij,ujij £ R such that 


( 2 ) 


Zj — (cr tj + iuJij)(zj Zj ). 


The scalars < 7,; 7 and describe the infinitesimal scalings and rotations of the edges. 
They satisfy the following compatibility conditions: 


Lemma 3.5. Given aij,u)ij £ R the following statements are equivalent: 

(a) There exist ii such that |I|) holds. 

(b) We have 

(3) 0 = (<T 12 + iu 12 )(z 2 - Zl) + (cr 23 + IU)23)(Z3 - z 2 ) + (031 + *w 3 i)(zi - Z3). 

(c) There exists u> £ R such that 

ioo = iuj 23 + i cot /3 i(ct 3 i - er 12 ) 

= iw 3 i + i cot /3 2 (cti 2 - cr 23 ) 

= ?Wi 2 + i cot /3 3 (cr 23 - cr 3 i). 

(d) There exist a £ R such that 

a = c 23 + i cot Pi (zw 3 i — iui 12 ) 

= 031 +i cot/3 2 («wi2 - iw 23 ) 

= CTi 2 + i cot p 3 (lU>23 - *w 3 i). 

Proof. The relation between (a) and (b) is obvious. We show the equivalence be¬ 
tween (b) and (c). With A denoting the signed triangle area we have the following 
identities: 


0 = (*0i - Zi),Zj - Zi), 

2A = (i(zj - Zi),z k - Zj), 

(i(zj - z^, i(zj - zf)) = (zj - Zi, Zj - Zi). 

Using these identities and Z3 — Z2 £ span R {?'(^i — z 3 ),i(z2 — Z\)} we obtain 
(4) z 3 - z 2 = cot(/ 3 3 )i(z 2 - Zi) - cot(/ 3 2 )i(zi - z 3 ). 

Cyclic permutation yields 

zi - z 3 = cot(/ 3 i )i(z 3 - z 2 ) - cot (/ 3 3 )i(z 2 - 21), 

Z2 - zi = cot(ft 2 )i(zi - z 3 ) - cot(/ 3 i )i(z 3 - z 2 ). 
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Substituting these identities into Equation (|3| we obtain 

0 = Oi ( C0t(/3 3 )f(z 2 - Zi) - COt(/? 2 )*(~l - z 3 )) + L 0 23 i ( z 3 ~ Z 2 ) 

+ er 2 ( cot(/3i)i(z 3 - z 2 ) - cot {/ d 3 ) i ( z 2 - Zi)) + w 3 i*(zi - z 3 ) 

+ cr 3 ( cot(/3 2 )i(zi z 3 ) -cot(/3i)*(z 3 -z 2 )) +wi 2 *(z 2 - zi) 

= ( u>i + cot / 3 i ( a 2 - < j 3 )) i ( z3 - z 2 ) 

+ (w 2 + cot /3 2 (ct 3 - t 7 i))i(zi - z 3 ) 

+ (w 3 + cot/3 3 (cri - cr 2 ))i(z 2 - zi). 

Now we use that Ai, A 2 , A 3 G C satisfy 

Ai*(z 3 - z 2 ) + A 2 ?'(zi - z 3 ) + A 3 i(z 2 - zi) = 0, 

if and only if Ai = A 2 = A 3 . This establishes the equivalence of (b) and (c). The 
equivalence of (b) and (d) is seen in a similar fashion by eliminating i(zj — zj) in 
([3]) instead of (Zj — zi). 

□ 

The quantity w above describes the average rotation speed of the triangle. Sim¬ 
ilarly, it can be verified that the above a satisfies 

R 



where R denotes the circumradius of the triangle. Thus er signifies an average 
scaling of the triangle. 

3.2. Harmonic functions with respect to the cotangent Laplacian. In smooth 
complex analysis conformal maps are closely related to harmonic functions. If a 
conformal map preserves orientation it is holomorphic and satisfies the Cauchy Rie- 
mann equations. In particular, its real part and the imaginary part are conjugate 
harmonic functions. Conversely, given a harmonic function on a simply connected 
surface then it is the real part of some conformal map. 

A similar relationship manifests between discrete harmonic functions (in the 
sense of the cotangent Laplacian) and infinitesimal deformations of triangular meshes. 
Discrete harmonic functions can be regarded as the real part of holomorphic func¬ 
tions which satisfies a discrete analogue of the Cauchy Riemann equations. In 
particular, a relation between discrete harmonic functions and infinitesimal pattern 
deformations was found by Bobenko, Mercat and Suris [2]. Integrable systems were 
involved in this context. We extend their result to include the case of infinitesimal 
conformal deformations. 

Theorem 3.6. Let z : V — > C be a simply connected triangular mesh realized in 
the complex plane and h : V —» R be a function. Then the following are equivalent: 

(a) h is a harmonic function with respect to the cotangent Laplacian, i.e. using 
the notation of Figure^ 7} for all interior vertices i £ we have 

(5) 5Z( cot $3 + cot )(hj ~ hi) =0. 

j 

(b) There exists an infinitesimal conformal deformation z : V — > C with scale 
factors given by h. It is unique up to infinitesimal rotations and translations. 

(c) There exists an infinitesimal pattern deformation iz : V —> C with h as 
angular velocities. It is unique up to infinitesimal scalings and translations. 
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Proof. We show the equivalence of the first two statements. The equivalence of the 
first and the third follows similarly. 

Suppose h is a harmonic function. Since the triangular mesh is simply connected, 
equation ©> implies the existence of a function Co : F —> R. such that for all interior 
edges {ij} we have 

'/.(Uj.y i,- lujjn — i (cot f3fj T cot {3j.f)(hj h t ). 


Here Co is unique up to an additive constant and called the conjugate harmonic 
function of h. Using Co we define a function oo : E —> R via 


UO>i 


uij — ZUJijk 7 COt 'd,'y (hj hi). 

Lemma 3.5 now implies that there exists z : V —> C such that 


(Zj - Zi) = 


hi -f- hj 


+ iijJij 


(Zj - Zi). 


This gives us the desired infinitesimal conformal deformation of z with h as scale 
factors. 

To show uniqueness, suppose z, z! are infinitesimal conformal deformations with 
the same scale factors. Then z — z' preserves all the edge lengths of the triangular 
mesh and hence is induced from an Euclidean transformation. 

Conversely, given an infinitesimal conformal deformation i with scale factors h. 
We write 

hi hi 


for some to : E 
that 

We have 
and 


+ iuijj (Zj - Zi) 

l. Lemma |3.5| implies that there is a function Co : F 
iiOjjk — ibOij T 7 cot Sij ( hj hi). 

iCoijk - iCoju = i(cot fiij + cot p l ji){hj - hi) 
£(cot4 + cotph)(hj - hi) = 0. 


such 


Therefore h is harmonic. 


□ 


4. HOLOMORPHIC QUADRATIC DIFFERENTIALS 

In this section, we introduce a discrete analogue of holomorphic quadratic dif¬ 
ferentials. We illustrate their correspondence to discrete harmonic functions. It 
reflects the property in the smooth theory that holomorphic quadratic differentials 
parametrize Mobius structures on Riemann surfaces (Ch. 9, ®)- 

To simplify the notation, we make use of discrete differential forms . We denote 
E the set of oriented edges and E int the set of oriented interior edges. Given an 
oriented triangular mesh M, a complex-valued function 77 : E —t C is called a 
discrete 1-form if 

hi^ij) = —r]{eji) We.ij G E. 

It is closed if for every face {ijk} 

v( e ij) + v( e jk) + r](eki) = 0 . 

It is exact if there exists a function / : V —> C such that 

vi^ij) df(eij) . fj — fi . 

Similarly, we can consider discrete 1-forms on the dual graph M* of M and these are 
called dual 1-forms. Given an oriented edge e, we denote e* its dual edge oriented 
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from the right face of e to its left face. The set of oriented dual edges is denoted 
by E*. 

Definition 4.1. Given a triangular mesh z : V C realized on the complex 
plane, a function q : Ej nt —> *K. defined on interior edges is a discrete holomorphic 
quadratic differential if it satisfies for every interior vertex i G V int 

/* ' Qij = 0 ) 

3 

'y Qij/dz(eij) = 0 . 

3 

Theorem 4.2. Let q : Ei nt zR. be a holomorphic quadratic differential on a 
realization z : V —> C of a triangular mesh. Suppose $ : C —> C is a Mobius trans¬ 
formation which does not map any vertex to infinity. Then q is again a holomorphic 
quadratic differential on w := 4> o 2 . 


Proof. Since Mobius transformations are generated by Euclidean transformations 
and inversions, it suffices to consider the inversion in the unit circle at the origin 

w := $( 2 ) = 1/z. 


We have 

y ' Qij/dw(eij) = y ' —ZiZjqij/dz(eij) = —Zi y ( qij — Zi y ( qij/dz{eij) = 0 . 

3 3 3 3 

Hence the claims follow. □ 


We are going to show that on a simply connected triangular mesh, there is a 
correspondence between discrete holomorphic quadratic differentials and discrete 
harmonic functions. 

We first show how to construct a discrete holomorphic quadratic differential from 
a harmonic function. Given a function u : F -> R on a realization of z : V —> C 
of a triangular mesh M. If we interpolate it piecewise-linearly over each triangular 
face, its gradient is constant on each face and we have grad, u : F —> C given by 


grad, Uijk = i 


Uidz{ejk) + Ujdz(e ki ) + u k dz(eij ) 


2 A. 


ijk 


Note that we ignore here the non-generic case (which leads to the vanishing of 
the area) where the triangle degenerates in the sense that its circumcircle passes 
through the point at infinity. Also note that for a non-degenerate triangle that is 
mapped by 3 in C in an orientation reversing fashion the area A l]k is considered to 
have a negative sign. Granted this, one can verify that the gradient of u satisfies 


(grad. u ijk , dz(eij )) = Uj - \/{ij} C {ijk} G F. 

We define u z : F —> C by 

u z ■■= 2§ rac b u - 

and the dual 1-form du z : E* nt C on M by 

du z (e i j) . (n z }ijk (u z )ju 

where {ijk} is the left face and {jil} is the right face of the oriented edge e z j. 


Lemma 4.3. Given a function u : V —> R on a realization of a triangular mesh 
z : V —> C, we have 


du z {e* ij )dz{e ij ) 

—X ■ . 

= —( cot/3 l jk (u k - Uj) + cot p ki (u k - Uj) + cot (3 3 u (ut - uf) + cot - uf)). 
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which is purely imaginary (Figure [7p. 

Proof. Since 

(gr&d z u ijk ,dz(eij)) = Uj - Ui = (grad z Uj k i,dz(eij)), 

we have 

R e(du z (e*j)dz(e.ij)) = 0. 

On the other hand, using equation Q we get 
R e(du z (e*j)idz(eij)) 

= R,e (((u 2 ) i-j k (u z )ju)idz(eij)) 

=((grad z Uijk, cot f3) k dz{e jk ) - cot f3 ki dz(e ki )) 

+ (grad, Uju, cot /3 J u dz{eu) - cotf}\ j dz(e lj )))/2 

= ^(cot/3j fc (u fe ~Uj)+ cot /3 J ki (u k - m) + cot /3^ (u t - uf) + cot /3y («; - uj)). 

Hence the claim follows. □ 

Lemma 4.4. Given a realization z : V —> C of a triangular mesh. A function 
u : V —> R is harmonic if and only if the function q : E int —> defined by 

Qij • du z (e i j)dz{eij) 
is a holomorphic quadratic differential. 

Proof. Note q is well defined since 

qij = du z (e* j )dz(e ij ) = du z (e* ji )dz(e ji ) = q jt . 

It holds for general functions u : V —> R. that 

Re(<?) = 0 

Y Qij/ dz ( e ij) = Y du z( e *j ) =° V* G V int . 

3 3 

We know from Lemma |4.3| that for every interior vertex i G Vi n t 

Y = du z {e*j)dz(e ZJ ) = ^(cot + cot P^ffuj - Ui ). 
j j j 

Hence, u is harmonic if and only if q is a holomorphic quadratic differential. □ 

Lemma 4.5. Let z : V —> C be a realization of a simply connected triangular mesh. 
Given a function q : Ei nt such that for every interior vertex i £ Vi n t 

Y<hj/ dz (^j) = 

3 

there exists a function u : V —> R such that for every interior edge {ij} 

qij = du z {e* ij )dz{e i j). 

Proof. We consider a dual 1-form r on M defined by 

r(e*j) = qij/dz{eij). 

Since M is simply connected and 

Y r ( e i?) = 5Z Qij/ dz ( e ij ) = °> 

3 3 

there exists a function h : F —> C such that 

dhfeij ) .— hijk hju — r(e^). 
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It implies we have Re(dh(e*)dz(e)) = Re(g) = 0 and 

•— (2/ijjfe, dz[e.ijf) — (fihjn, dz^e-jj)). 

is a well-defined R-valued 1-form. Since the triangular mesh is simply connected 
and for every face {ijk} 

) “t - aj(ej k ) H- aj(e k f) — 0, 

there exists a function u : V -* R such that for every oriented edge e^- 

du(e.ij) = Uj — Ui= u(eij). 

It can be verified that 

h = -grad z M = u z . 

Hence we obtain 

Qij = T ( e i j )dz(eij) = dh(e* j )dz(ei j ) = du z {e* ij )dz{e i j) 
for every interior edge {ij}. □ 


Theorem 4.6. Suppose z : V —> C is a realization of a simply connected triangular 
mesh. Then any holomorphic quadratic differential q : E int —> iR is of the form 

—V 

qij — du z (eij)dz(eij) if (E E ?rl f 
for some harmonic function u : V —> M. 

Furthermore, the space of holomorphic quadratic differentials is a vector space 
isomorphic to the space of discrete harmonic functions module linear functions. 


Proof. The first part of the statement follows from Lemma 4.4 and Lemma 4.5 In 
order to show the second part, it suffices to observe that 


du z = 0 
for some a, b £ C. 


grad u = a 


du = (a, dz) 


u = (a. z 


)+b 


□ 


In previous sections, we showed that every harmonic function corresponds to an 
infinitesimal conformal deformation. The following shows that discrete holomorphic 
quadratic differentials are the change in the intersection angles of circumscribed 
circles. 


Theorem 4.7. Let z : V —> C be a realization of a simply connected triangular 
mesh. Suppose u : V —> M is a discrete harmonic function and z is an infinitesimal 
conformal deformation with u as scale factors. Then we have 


du z dz = — 


1 c'r 2 

2 cr. 



where : E int —>• R denotes the change in the intersection angles of neighboring 
circumscribed circles. 


Proof. We write (ij — if) = ( h, + h i -g iunj)(zj — zf). Applying Lemma 
ci' z jj l cr z jj z izjj i,. icoki 4~ ioJu iujij 


4.3 


we have 


=i(cot/3j k (u k - Uj) + cot/3 J ki (u k - uf) + cot - uf) + cot l3\j(ui - Uj )) 
= - 2 du z {e* ij )dz(e i j). 


The equality 



follows from Equation ([T]). □ 
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5. Conformal deformations in terms of End(C 2 ) 

In this section we show how an infinitesimal conformal deformation gives rise to 
a discrete analogue of a holomorphic null curve in C 3 . Later we will see that the 
real parts of such a ’’holomorphic null curve” can be regarded as the Weierstrass 
representation of a discrete minimal surface. 

Up to now we have mostly treated the Riemann sphere CP 1 as the extended 
complex plane C = C U {00} . In this section we will take a more explicitly Mobius 
geometric approach: We will represent fractional linear transformations of C by 
linear transformations of C 2 with determinant one. Actually, the group of Mobius 
transformations is 

Mob(C) “ PSL(2,C) “ SL(2,C)/(±J). 

However, since we are mainly interested in infinitesimal deformations and any map 
into PSL(2,C) whose values stays close to the identity admits a canonical lift to 
SL(2,C), we can safely ignore the difference between PSL(2,C) and SL(2,C). 

Given a realization z : V —> C of a triangular mesh we consider its lift if : V —> C 2 



and regard the realization as a map T : V —> CP 1 where 

* := C ( l ) = [# 

The action of a Mobius transformation on the Riemann sphere is given by a 
matrix A £ SL(2, C), which is unique up to sign: 

M ^ [Av\- 

Before we investigate infinitesimal deformations we first consider finite deforma¬ 
tions of a triangular mesh : V —> CP 1 . Given such a finite deformation, the 
change in the positions of the three vertices of a triangle {ijk} can be described by 
a Mobius transformation, which is represented by Gijk £ SL(2, C). They satisfy a 
compatibility condition on each interior edge {ij} (see Figure [I]): 

[Gijkipi] = [Gjuipi], 

[Gijklpj] = [Gjillpj]- 

Suppose now that the mesh is simply connected. Then up to a global Mobius 
transformation the map G : F —> SL(2,C) can be uniquely reconstructed from the 
multiplicative dual 1-form defined as 

G(eT) := G y \G ijk . 

G(e*j) is defined whenever {ij} is an interior edge and we have 

G(e*j) = G(e* i )~ 1 . 

Moreover, for every interior vertex i we have 

II G(e* j ) = I. 

j 

The compatibility conditions imply that for interior each edge {ij} there exist 
£ C\{0} such that 

G{e*j)ipi = Xij'iipi 

= 'Vj. 
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Since A ijj A ijj = det(G(e* ) )) = 1, we have 

A ij • \ij t i 1 /Xyj.j ■ 

Because of G(e*j) = G(e* J ) _1 we know 

A ij Xjjj 1/A jii A ji~ 


Hence A defines a complex-valued function on the set Ei nt of interior edges. 

We now show that for each interior edge A ^ determines the change in the cross 
ratio of the four points of the two adjacent triangles. Note that the cross ratio of 
four points in CP 1 can expressed as 


cr{[ipj\, [ipi c ], [ipi], [ipi]) 


det(ip k ,ipj)det{ipi,ipj) 

det(ipi,ip k )det(ipj,ipi)' 


Lemma 5.1. Suppose we are given four points [ipi], [ ipj ], [ip k ], [ipi] G CP 1 and G £ 
SL(2,C) with 

Gifi = A 


Gipj = Xipj 


for some A £ C\{0}. Then the cross ratio of the four transformed points 

W>i] = [Gtpi] , [ipj] = [Gipj] , [ipk] = [Gip k ] , [ipi] = [ipi] 

is given by 


Proof. 


cr ([Vh']> bPk], bPi], bPi]) = cr ([ipj], [ip k \, bPi], bPi])/x 2 . 


cr ($j], [V’fc], bPi], [in]) 


det (Gip k , Gipj) det (ipi, Gipj) 
det ( Gifi, Gip k ) det ( Gipj, ipi ) 
cr ([ipj], [ip k ], [ipi ], [ipi])/X 2 . 


□ 


We now can summarize the information about finite deformations of a realization 
as follows: 


Theorem 5.2. Let d* : V — > CP 1 be a realization of a simply connected triangular 
mesh. Then there is a bijection between finite defoliations of d* in CP 1 modulo 
global Mobius transformations and multiplicative dual 1 forms G : E* nt -* SL(2,C) 
satisfying for every interior vertex i 

n G <4)= 7 

j 

and for every interior edge 

G{e*j) = G(e*,) _1 
G(e*j)iPi = Ar-Vi 
G(e*j)ipj = Xijipj. 

Here X : Ei nt — > C\{0}. We denote by cr : Ei nt C the cross ratios of and 
cr : E int — > C the cross ratios of a new realization described by G. Then 

cr = cr / A 2 . 


In particular, 


| A | = 1 => the deformation is conformal. 

Arg(A) = 0 => the deformation is a pattern deformation. 



14 


WAI YEUNG LAM AND ULRICH PINKALL 


Suppose we have a family of deformations described by dual 1-forms G* : Ei nt —> 
SL(2, C) with Go = I. By considering rj := pjp\t=o G t we obtain the following 
description of infinitesimal deformations: 

Corollary 5.3. Let d* : V —> CP 1 be a realization of a simply connected triangular 
mesh. Then there is a bijection between infinitesimal deformations of d/ in CP 1 
modulo infinitesimal Mobius transformations and dual 1 forms p : Ei nt -A sl(2, C) 
satisfying for every interior vertex i 

( 6 ) Y^r)(et j ) = 0 

3 

and for every interior edge 

nie-ij) = -v(e*ji ) 

V(e*j)i>i = ~Tij i>i 
v(e*j)ipj = IMj ipj- 

Here /i : Ei nt —> C. We denote by cr : Ei nt -A C the cross ratios of d/ and cr : 
Ei n t -A C the rate of change in cross ratios induced by the infinitesimal deformation 
described by rj. Then 

1 cr 
» = ~2n 

In particular, 

Re(/r) = 0 ==> the infinitesimal deformation is conformal , 

Im(/i) = 0 => the infinitesimal deformation is a pattern deformation. 

Note that given a mesh, the 1 -form ij is uniquely determined by the eigenfunction 
p. We now investigate the constraints on p implied by the closedness condition Q 
of 77 . 

Consider the symmetric bilinear form (,):C 2 xC 2 A sl(2,C) 

(<j>, <p)v := det(<p, v)(p + det(y>, v)<j>. 

For ipi 7 ^ ipj £ C 2 we define 

The matrix mij is independent of the representatives of [ipi], [ipj] £ CP 1 and we 
have 


Tflij — TTlj i 

niijipi = ipi 

mijipj = ipj. 


Using the representatives ipi = 
V( e ij) = 


Zi 

1 

Pij 


det {ipj, ipi 


we obtain 


-(ipi,ipj 


Tij ( Zi + Zj 
Zj - Zi\ 2 


2 Zi Zj \ 

-U: - Zj J 


Hence 


( 7 ) 


Y ^* 3 ) = 0 <*=>• Y Vij = 0 and Y Vij / ( z j - Zi) = 0. 

3 3 3 
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We consider the Pauli matrices 



which form a basis of sl(2,C). Then 

= J 7 ((1 - ZiZj)ai + i(l + ZiZj)a 2 + (zi + zj)a 3 ). 

Zj Zi 

If we now identify sl(2, C) with C 3 via 

/ 1 \ 

tJi 0 , cr 2 

V 0 ) 

we obtain 


1 - ZiZj \ 

i{l + ZiZj) 

Zi + Zj ) 




Thus to every infinitesimal deformation of a realized triangular mesh we can as¬ 
sociate a closed sl(2, C)-valued dual 1-fornr. In the special case of an infinitesimal 
conformal deformation (i.e. /r is real-valued) we will see that this yields a discrete 
analogue of the Weierstrass representation for minimal surfaces. 


6. Weierstrass representation of discrete minimal surfaces 


The Weierstrass representation for minimal surfaces in R 3 is the most classical 
example for applications of complex analysis: 

Theorem 6.1. Given two meromorphic functions g,h : U cC-)C such that g 2 h 
is holomorphic. Then /:£/—>■ R 3 defined by 

i-g 2 \ \ (( i-3 2 

*(! + g 2 ) I h(z)dz\ = Re I I i(l + g 2 ) 

is a minimal surface. Its Gaufl map n is the stereographic projection of g 

2 Re g \ 

2 Im g I . 

Isl 2 -i / 

The holomorphic quadratic differential q := hg z dz 2 is called the Hopf differential 
of f and encodes its second fundamental form: The direction defined by a nonzero 
tangent vector W is 

an asymptotic direction <=> qiW) £ ?'R. 
a principal curvature direction q(W) £ R. 

Locally, every minimal surface can be written in this form. 




df = Re 


We now develop a discrete version of this theorem for arbitrary triangular meshes 
realized in the complex plane. A similar formula for quadrilateral meshes with 
factorized real cross ratios was established by Bobenko and Pinkall [T]. Here we 
will use the definition of a discrete minimal surface / with Gaufi map n given in 
0 : 

Definition 6.2. Let n : V —> S 2 be a realization of a triangular mesh on the unit 
sphere in R 3 . Then a map f : F R 3 defined on the set F of faces is called a 
discrete minimal surface with GauB map n if for all oriented interior edges e l3 we 
have 

(jlj II?!) X (fijk fjil ) 0’ 





16 


WAI YEUNG LAM AND ULRICH PINKALL 



Figure 3. Left: a triangulated surface n : V —► § 2 with vertices 
on the unit sphere. Right: a discrete minimal surface / : F —> R 3 
satisfying Definition |6.2| 


Here {ijk} and {jil} denote the left and the right faces of e,j. 

This definition mirrors the fact from the smooth theory that minimal surfaces 
are Christoffel duals of their Gaufi maps (Figure [3]). The correspondence between 
discrete harmonic functions and discrete minimal surfaces was observed in [5]. Here 
is a Weierstrass representation for discrete minimal surfaces in terms of their GauB 
map and their Hopf differential: 


Theorem 6.3. Let z : V —► C be a realization of a simply connected triangular 
mesh and q : Ei nt A # a holomorphic quadratic differential. Then there exists 
f : F —> R 3 such that for every interior edge {ij} 


(9) 


df{etj) = Re 



1 - ZiZj \ \ 
*(1 + Z,Zj) 

Zi +Zj j j 


Moreover f is a discrete minimal surface with Gaufi map 


n = 


1 

l?Ti 


2 R ez \ 
21 Ill 2 I 

l*| 2 -l J 


Locally, every discrete minimal surface can be written in this form. 


Proof. Suppose q : E int — >• zM is a holomorphic quadratic differential. Then by 0 
and Q the dual 1 -form 77 defined as 


V(e*j) ■= 


Qij 

i{zj - zf) 


1 - ZiZj \ 
i(l + ZiZj) 

Zi + Zj J 


satisfies 

=0 ' 

j 

for all interior vertices i. Therefore, since the triangular mesh is simply connected, 
there exists 3 : F —> C 3 such that for any interior edge e we have 


d${e*) = V {e*). 

Thus the map / : F —► M 3 defined by f := Red$ satisfies Equation 0 - To show 
that / is a discrete minimal surface we define a function k : E lnt — > R by 

kij ■.= —i qij/\zj - Zi\ 2 . 
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Then by direct computation we obtain 


( 10 ) 


df(e*A = 


Mi + N 2 )(i + M 2 ) 


(nj — rii). 


~ 2 v 3 

This shows that / is a discrete minimal surface with Gaufi map n. The converse 
is straightforward: Given a discrete minimal surface / with GauB map n we define 
k : Ei n t R via (10). Then it can be shown that the function 

Qij • — 7 k-'ij | Zj Zi | 

is a holomorphic quadratic differential. □ 


Remark 6.4. The discrete minimal surfaces given by ([9]) are trivalent meshes 
with planar vertex stars for purely imaginary q. It is closely related to discrete 
asymptotic nets. The factor i in front of z 3 — Zi appears since the integration is 
taken over a dual mesh while in the smooth theory *dz = idz. 


Note that we could also consider the periodic one-parameter family of maps 
/“ : F —> R 3 defined for a £ R by 

r := Re(e iQ 5). 

This family of discrete surfaces can be regarded as an associate family of minimal 
surfaces and is investigated in J7j. 


References 

[1] A. Bobenko and U. Pinkall. Discrete isothermic surfaces. J. Reine Angew. Math., 475:187- 
208, 1996. 

[2] A. I. Bobenko, C. Mercat, and Y. B. Suris. Linear and nonlinear theories of discrete analytic 
functions. Integrable structure and isomonodromic Green’s function. J. Reine Angew. Math., 
583:117-161, 2005. 

[3] A. I. Bobenko, U. Pinkall, and B. A. Springborn. Discrete conformal maps and ideal hyper¬ 
bolic polyhedra. Geom. Topol., 19(4):2155-2215, 2015. 

[4] R. J. Duffin. Basic properties of discrete analytic functions. Duke Math. J., 23:335-363, 1956. 

[5] J. Ferrand. Fonctions preharmoniques et fonctions preholomorphes. Bull. Sci. Math. (2), 
68:152-180, 1944. 

[6] R. C. Gunning. Lectures on Riemann surfaces. Princeton Mathematical Notes. Princeton 
University Press, Princeton, N.J., 1966. 

[7] W. Y. Lam. Discrete minimal surfaces: critical points of the area functional from integrable 
systems. 2015. arXiv: 1510.08788. 

[8] W. Y. Lam and U. Pinkall. Isothermic triangulated surfaces. 2015. arXiv: 1501.02587. 

[9] F. Luo. Combinatorial Yamabe flow on surfaces. Commun. Contemp. Math., 6(5):765-780, 
2004. 

[10] C. Mercat. Discrete Riemann surfaces and the Ising model. Comm. Math. Phys., 218(1): 177— 
216, 2001. 

[11] O. Schramm. Circle patterns with the combinatorics of the square grid. Duke Math. J., 
86(2):347-389, 1997. 

[12] S. Smirnov. Discrete complex analysis and probability. In Proceedings of the International 
Congress of Mathematicians. Volume I, pages 595—621. Hindustan Book Agency, New Delhi, 
2010 . 

[13] B. Springborn, P. Schroder, and U. Pinkall. Conformal equivalence of triangle meshes. ACM 
Transactions on Graphics (TOG), 27(3):77, 2008. 

Wai Yeung Lam, Technische Universitat Berlin, Institut fur Mathematik, Strasse 
des 17. Juni 136, 10623 Berlin, Germany 

Ulrich Pinkall, Technische Universitat Berlin, Institut fur Mathematik, Strasse 
des 17. Juni 136, 10623 Berlin, Germany 

E-mail address: lam@math.tu-berlin.de, pinkall@math.tu-berlin.de 




